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ABSTRACT 

We present time resolved J-band (1.025 - 1.340 fim) infrared spectra of the short- 
period dwarf novae (DNe) WZ Sge, VY Aqr and single spectra of the short-period 
DNe EF Peg and the novalike variable PX And. There is some evidence in the spectra 
of VY Aqr and EF Peg that we have detected the secondary star, both in the con- 
tinuum slope and also through the possible presence of spectral features. The spectra 
of WZ Sge and PX And, on the other hand, show no evidence for the secondary 
star, with upper limits for it's contribution to the J-band light of 10% and 20% re- 
spectively. The spectral type of the secondary in WZ Sge is constrained to be later 
than M7.5V. Using skew mapping we have been able to derive a value for the radial 
velocity semi-amplitude of the secondary star in VY Aqr of Kr — 320 ± 70 kms -1 , 
which in conjunction with K\y from Thorstensen & Taylor (1997) gives a mass ratio 
of q = 0.15 ±0.04. 

Key words: binaries: close - stars: individual: WZ Sge, VY Aqr, EF Peg, PX And 
- dwarf novae, cataclysmic variables - infrared: stars 



1 INTRODUCTION 

Cataclysmic variables (CVs) are semi-detached binary sys- 
tems in which a white dwarf primary accretes material from 
a Roche-lobe filling secondary. For a thorough review of 
CVs see Warner (1995). There is a strong motivation to 
detect the secondary stars in CVs. A detection of the sec- 
ondary star is desirable for a reliable estimation of the sys- 
tem parameters, which in turn allow one to test theories 
of CV evolution, (e. g. see Smith & Dhillon 1998). In ad- 
dition the contribution of the secondary to the total flux 
can be used to give the distance to the CV (e. g. see Marsh 
1990). Detection of the secondary star is of particular in- 
terest in CVs below the period gap. The models of Kolb 
(1993) predict that 99% of the present-day population of 
CVs should be below the period gap and that around 70% 
of these systems will have already reached the orbital pe- 
riod minimum. Modelling of the secondary stars in post- 
period-minimum CVs has been performed by Howell, Rap- 
paport & Politano (1997), who find them to be degener- 
ate, brown dwarf-like objects with masses between 0.02-0.06 
Mq. Howell, Rappaport & Politano (1997) further specu- 
late that the tremendous outburst amplitude dwarf novae 
(TOADs) are these post-period-minimum CVs. Estimates 



of the secondary stars' parameters would provide us with 
tests for these claims, as well as allowing us to investigate 
the mass-radius relation of the brown dwarf-like secondary 
by invoking the Roche lobe-filling criterion. 

Previous attempts to detect the secondary star in CVs 
have focussed on the TiO bands and Nal lines in the optical- 
infrared (0.7 - 1.0 pm), e.g Friend et al. (1988). These 
searches proved insensitive to the secondary star in short- 
period CVs, although there were exceptions (e.g Z Cha, 
Wade & Home 1988). The relative insensitivity of optical 
searches to the secondary in short-period CVs is unsurpris- 
ing, as short period CVs are expected to possess late-type 
M-dwarf secondary stars*. These late- type M-dwarf stars 
have spectra which peak in the infrared (see Kleinmann & 
Hall 1986), and at optical wavelengths the light from the 
M-dwarf is swamped by light from the accretion disc. In- 
frared searches for the secondary star in short period CVs 
have largely focussed on the K-band (Dhillon 1997), with 



* The spectral type of the secondary star in a CV can be 
estimated from it's orbital period, P, using the relationships 
26.5 - 0.7P (for P < 4 hr) and 33.2 - 2.5P (for P > 4 hr), 
where GO = 0, K = 10 and M0 = 20 (Smith & Dhillon 1998) 
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some success. Whilst the secondary star has proved rela- 
tively simple to find in dwarf novae above the period gap, 
where the contribution of the secondary in the K-band is 
typically > 75% (Dhillon 1997), it has proved much more 
difficult to detect the secondary star in both nova-likes and 
dwarf novae below the period gap. For these systems the 
contribution of the secondary to the infrared flux is usually 
much less (typically less than 25%; Dhillon 1997). To find 
the secondary stars in these systems a new strategy is re- 
quired. In this paper we present the first results of a program 
to detect secondary star features using the J-band spectra 
of cataclysmic variables. 

Searching for secondary star features in the J-band has 
certain advantages over other wavelength ranges. First, ob- 
servations in the near-infrared are plagued by a very high 
background. This is due to both thermal emission from the 
telescope, sky and instrument, and atmospheric emission 
from molecules (mostly OH~ and O2) excited by daytime 
solar radiation, resulting in airglow which is both temporally 
and spatially variable (Ramsay, Mountain & Geballe 1992). 
This background is significantly lower in the J-band than 
in K (16 mag/arcsec in J compared with 12 mag/arcsec in 
K - data for Mauna Kea, Hawaii). Second, in the late-type 
M-dwarfs expected in post-period gap CVs, there are more 
absorption features in the J-band than K, and these absorp- 
tion features have larger equivalent widths than those in the 
K-band (e.g see Jones et al. 1994). Third, the change in 
strength of absorption features with spectral type in late- 
type M-dwarfs is more marked in the J-band than in K (c.f. 
table 2 in Dhillon & Marsh 1995 with table 2 in this paper). 
This high sensitivity of absorption features to spectral type 
makes a secondary star detection in the J-band a power- 
ful constraint on models of the evolution of CVs below the 
period gap. 



2 OBSERVATIONS 

On the nights of 1998 August 8 & 9 we obtained spectra 
of the Dwarf Novae (DNe) WZ Sge, EF Peg, VY Aqr, the 
novalike variable PX And and the M-dwarfs G1644C and 
G1866AB with the Cooled Grating Spectrometer (CGS4) on 
the 3.8 m United Kingdom Infrared Telescope (UKIRT) on 
Mauna Kea, Hawaii. CGS4 is a 1-5 micron spectrometer 
containing an InSb array with 256x256 pixels. The 40 1/mm 
grating with the 300 mm camera gave a resolution of approx- 
imately 300 kms" 1 . To cover the wavelength range 1.025- 
1.340 microns required one grating setting, centred at 1.175 
microns (second order) . Optimum spectral sampling and bad 
pixel removal were obtained by mechanically shifting the ar- 
ray over two pixels in the dispersion direction in steps of 0.5 
pixels. We employed the non-destructive readout mode of 
the detector to reduce the readout noise. In order to com- 
pensate for fluctuating atmospheric OH~ emission lines we 
took relatively short exposures (typically 60s) and nodded 
the telescope primary so that the object spectrum switched 
between two different spatial positions on the detector. The 
slit width was 0.6 arcseconds (projecting to approximately 
1 pixel on the detector) and was oriented at the parallactic 
angle throughout the second night. Observations of WZ Sge 
and G1866AB were taken with a slit position angle of 0° to 
avoid light from a star close to WZ Sge falling on the slit. 



Using the UKIRT tip-tilt secondary, the seeing was around 
0.5" on the first night, and dropped steadily from 0.8" to 
0.4" on the second night. The humidity was low (10-20%) 
throughout, and the sky was photometric for most of the 
run, although some high cirrus developed in the second half 
of the first night. A full journal of observations is presented 
in table 1. 



3 DATA REDUCTION 

The initial steps in the reduction of the 2D frames were per- 
formed automatically by the CGS4 data reduction system 
(Daley & Beard 1994). These were: the application of the 
bad pixel mask, bias and dark frame subtraction, flat field 
division, interlacing integrations taken at different detector 
positions, and co-adding and subtracting nodded frames. 
Further details of the above procedures may be found in 
the review by Joyce (1992). In order to obtain ID data, we 
removed the residual sky by subtracting a polynomial fit and 
then extracted the spectra using an optimal extraction tech- 
nique (Home 1986). The next step was the removal of the 
ripple arising from variations in the star brightness between 
integrations (i.e. at different detector positions). These vari- 
ations were due to changes in the seeing, sky transparency 
and the slight motion of the stellar image relative to the slit. 

There were two stages to the calibration of the spectra. 
The first was the calibration of the wavelength scale using 
krypton arc lamp exposures. The second-order polynomial 
fits to the arc lines yielded an error of less than 0.0001 mi- 
crons (rms). The final step in the spectral calibration was 
the removal of telluric features and flux calibration. This was 
performed by dividing the spectra to be calibrated by the 
spectrum of an F-type standard, with its prominent stellar 
features interpolated across. F-types were taken through- 
out the night, at different airmasses. In each case, the F- 
star used was that which gave the best removal of telluric 
features, judged by the cancellation of the strong features 
around 1.14 microns. We then multiplied the result by the 
known flux of the standard at each wavelength, determined 
using a black body function set to the same effective tem- 
perature and flux as the standard. As well as providing flux 
calibrated spectra, this procedure also removed telluric ab- 
sorption features from the object spectra. 



4 RESULTS 

Figure 1 shows the 1.025-1.340 /jm averaged spectra of the 
short period DNe WZ Sge and VY Aqr (left hand side) and 
the 1.025-1.340 fim averaged spectra of the short period 
DN EF Peg and the novalike variable PX And (right hand 
side) together with the spectra of M7 and M5 field dwarfs. 
We also show the spectrum of an F8V star, which indicates 
the location of telluric absorption features; spectral features 
within the strongest absorption bands are highly uncertain. 
The location of bad pixels are labelled by a bar showing the 
extent of the anomaly and the label 'bp'. In table 2 we list 
the wavelengths, equivalent widths and velocity widths of 
the most prominent spectral lines identified in figure 1. 
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Figure 1. Left-hand side: Average J-band spectra of the short period DNe WZ Sge, VY Aqr and the M-dwarfs G1644C and G1866AB. 
Right-hand side: Average J-band spectra of the short period DN EF Peg, the novalike variable PX And and the M-dwarfs G1644C and 
G1866AB. For both plots, the spectra have been normalized by dividing by the flux at 1.3 /im and then offset by adding a multiple of 
0.9 to each spectrum. Also shown is the spectrum of an F8V star, normalized by dividing by the flux at 1.3 (im, which indicates the 
location of telluric features. 



4.1 Accretion disc features 

The spectra of WZ Sge and VY Aqr (fig. 1) are dominated 
by the strong, double-peaked emission lines of Hel, Paschen- 
7 and Paschen-/3. Also seen is a small broad peak around 
1.07 /xm. This is identified as the Hel line at 1.0668 /im. 
The broad, double-peaked nature of these lines indicates an 
origin in the accretion disc. The full-width at half-maxima 
(FWHM) for WZ Sge are in agreement with those presented 
by Dhillon et al. (1999), who find a FWHM for the Brackett- 
7 line of 2200 kms" 1 . FWHM for VY Aqr are in good 
agreement with Thorstensen & Taylor (1997), who find a 
FWHM of 1550 kms -1 for Ha. We measured the separation 
of the double-peaked emission lines for VY Aqr and WZ Sge 
For WZ Sge the peak separations were 1400 ± 100 kms -1 
1300±100 kms" 1 and 1400±100 kms" 1 for Hel, P7 and P/3 
respectively. These values are in good agreement with val- 
ues for the peak separation derived from optical data, e.g 
1380 kms" 1 for Ha (Gilliland, Kemper & Suntzeff 1986) 
For VY Aqr the peak separations were 700 ± 100 kms" 1 
800 ± 150 kms" 1 and 900 ± 100 kms" 1 for Hel, P7 and 
P/3, respectively. These peak separations for VY Aqr are in 



good agreement with Thorstensen & Taylor (1997), who find 
a peak separation of 930 kms" 1 for Ha. 

The spectrum of PX And (fig. 1) is dominated by the 
strong, single-peaked emission lines of Hel, Paschen-7 and 
Paschen-/3, typical of the SW Sex stars of which PX And 
is a member (Still, Dhillon & Jones 1995). The FWHM for 
PX And are consistent with the values given by Still, Dhillon 
& Jones (1995), who find FWHM of 930 - 1030 kms" 1 for 
Hel emission and 1080 - 1600 kms" 1 for the Balmer lines. 

The spectrum of EF Peg (fig. 1) also shows strong emis- 
sion lines, although the signal-to-noise ratio is too low to 
say if these are single or double-peaked. The large velocity 
widths of the emission lines in EF Peg, suggests an accretion 
disc origin. 



4.2 Secondary star features 

The spectrum of VY Aqr shows increasing flux towards 
the red end of the spectrum, with a continuum slope of 
(7.4 ± 0.4) x 10" 1 mjy//tm, suggesting that the secondary 
star makes a significant contribution to the continuum light 
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Figure 2. Average J-band spectra of the short period DNc 
VY Aqr, EF Peg and the M-dwarf G1644C. The spectra have 
been normalized by dividing by the flux at 1.3 fim and then off- 
set by adding a multiple of 0.9 to each spectrum. One spectrum 
of G1644C has been divided by 1.7, in an attempt to match the 
appearance of the absorption feature at 1.33 (im to that seen in 
VY Aqr and noise has then been added to the spectrum in order 
to match the signal-to-noise to that of the DNe spectrum. Also 
shown is the spectrum of an F8V star, normalized by dividing 
by the flux at 1.3 fim, which indicates the location of telluric 
features. 

in VY Aqr (c.f. the slope of the M-dwarfs in figure 1). The 
continuum in VY Aqr also exhibits a change in slope at ~ 
1.1/Ltm. This feature is also present in the spectra of the 
M-dwarfs. The spectrum of VY Aqr also shows a tenta- 
tive detection of secondary star absorption features, par- 
ticularly the headless water band seen at 1.33 /im. Likewise, 
the EF Peg continuum shows increased flux at the red end of 
the spectrum, with a continuum slope of (5.7 ± 0.8) x 10 _1 
mjy//im and like the spectrum of VY Aqr, EF Peg also 
shows a tentative secondary star detection through the wa- 
ter band at 1.33 fj,m. The EF Peg continuum also shows a 
change in slope at ~ 1.1pm, in common with VY Aqr and the 
M-dwarfs. The detection of the water band in EF Peg and 
VY Aqr should be considered carefully: the spectra do not 
show any of the other absorption features visible in the spec- 
tra of the M-dwarfs - although it is possible they are lost in 



the noise. Figure 2 shows the spectra of the DNe EF Peg and 
VY Aqr, together with two spectra of the M-dwarf Gl 644c. 
One spectrum of Gl 644c has been divided by 1.7 in order 
to weaken the appearance of the spectral lines (the correct 
amount was judged by attempting to match the appearance 
of the 1.33/itm water feature to that seen in VY Aqr). This 
spectrum has then had noise added in order to match the 
signal-to-noise to that of the DNe spectra. Absorption fea- 
tures are barely discernable in this spectrum. This shows 
that absorption features from the secondary star are prob- 
ably hidden by noise in the CV spectra. Furthermore, the 
most prominent absorption features fall at the same wave- 
lengths as the strongest telluric features and bad pixels, and 
any features in these regions are inevitably more uncertain. 
In addition, the water feature at 1.33 /im only occurs in 
concurrence with the red continuum seen in VY Aqr and 
EF Peg, and does not appear in the spectrum of WZ Sge, 
for example. We believe that this constitutes a body of ev- 
idence for the detection of the secondary star in VY Aqr 
and EF Peg. This is discussed further in section 5.1. Un- 
fortunately, without detection of more absorption features 
and the observation of more field-dwarf templates we are 
unable to determine the spectral type of the secondary star 
in EF Peg and VY Aqr. 

There is no evidence for the secondary star in the spec- 
tra of WZ Sge and PX And, which show a blue continuum 
without any changes in slope or absorption features. 

4.3 Secondary star contributions and distances 

Table 3 shows the contribution of the secondary star to the 
total flux of each CV. Percentages were calculated by nor- 
malising the average spectra of both the CV and M-dwarf 
template. An increasing fraction of the M-dwarf spectrum 
was then subtracted from the CV spectrum until spectral 
features in the CV spectrum were removed or absorption fea- 
tures from the M-dwarf appeared in emission. The fraction 
at which this occurs gives the contribution of the secondary 
star to the total flux. Where no secondary star features are 
present this represents an upper limit to the contribution of 
the secondary star. 

The contributions were found to vary depending upon 
which M-dwarf template was used. It follows that systematic 
error could be introduced through using a template whose 
spectral type differs greatly from the actual spectral type of 
the secondary star. 

The strongest features in the template spectra increase 
in strength with spectral type. Hence templates of earlier 
spectral type give a higher contribution than those later 
in the spectral sequence. Hence, for WZ Sge and PX And, 
where the earliest spectral type consistent with the spectral 
type-orbital period relation in Smith & Dhillon (1998) was 
used, we have derived strict upper limits to the contribution 
of the secondary star to the J-band light. 

In VY Aqr and EF Peg, we used the M-dwarf which was 
able to correct the H2O band at 1.33 fim and flatten the con- 
tinuum without introducing spectral features in emission. It 
would have been apparent if the spectral type of the tem- 
plate used was too late, as the KI feature at ~ 1.17 /jm 
increases in strength more rapidly with spectral type than 
the H2O feature at 1.33 pm. Hence a template with too late a 
spectral type would not be able to correct the water feature 
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without introducing the KI feature in emission. Templates 
whose spectral types are too early would give a larger value 
for the secondary star contribution. However, the spectral 
type of the secondary in EF Peg and VY Aqr should be 
later than M5V. Hence we believe the contributions derived 
for these stars are reliable. The errors shown in table 3 are 
the discrepancies between the values obtained with the two 
M-dwarf templates and hence reflect the uncertainty result- 
ing from the unknown spectral type. 

The contributions of the secondary star will also be af- 
fected by systematic error due to flux from the disc if an 
emission feature corrupts the water band at 1.33 /tm. This 
is not thought to be the case. 

We find an upper limit to the secondary star contribu- 
tion to the J-band light in WZ Sge of 10%, consistent with 
Ciardi et al. (1998), who model the relative contributions 
to the near-infrared emission in WZ Sge and find that the 
secondary star should contribute w 20% of the near-infrared 
flux, although their models suggest that this figure is lower 
(perhaps « 10%) in the J-band. 

If the J-band magnitude for the CV is known it is pos- 
sible to derive the apparent J-band magnitude of the sec- 
ondary. In principle, this apparent magnitude can then be 
used to estimate the distance to the CV, by comparison 
with the absolute magnitude of a field-dwarf of appropri- 
ate spectral type. In the cases where the secondary has not 
been detected, this value represents a lower limit to the dis- 
tance involved. J-band magnitudes for WZ Sge, VY Aqr and 
PX And are 14.2, 15.3 and 14.2 respectively (Mark Huber, 
private correspondence) . The J-band magnitude of EF Peg 
is 16.1 (Sproats, Howell & Mason 1996). This gives a J- 
band magnitude of 15.7 ± 0.5 for the secondary in VY Aqr, 
16.7 ± 0.5 for the secondary in EF Peg and upper limits of 
16.2 and 15.4 for the J-band magnitudes of the secondaries 
in WZ Sge and PX And respectively. Unfortunately, we are 
unable to give distances to these CVs, as the absolute J-band 
magnitude of M-dwarfs varies rapidly with spectral type - 
c.f. Mj = 6.56 for M5V with Mj = 8.75 for M7V (Bessell 
1991) - a reliable spectral type determination is necessary 
before distances can be estimated. As an aside, it is worth 
noting that both Spruit & Rutten (1998) and Smak (1993) 
find a distance to WZ Sge of 48 pc, which, through the anal- 
ysis above, and the data in Bessell (1991) would constrain 
the spectral type of the secondary to be later than M7.5V, 
consistent with Ciardi et al. (1998), who find that the sec- 
ondary star in WZ Sge is cooler than 1700K. 

4.4 Time resolved spectra 

Figure 3 shows the J-band time-resolved spectra of VY Aqr 
(left-hand side) and the J-band time resolved spectra of 
WZ Sge (right-hand side) together with the spectra of 
M7 and M5 field dwarfs. The majority of ephemerides for 
WZ Sge are for mid-eclipse of what is believed to be the 
bright spot. In this paper, WZ Sge is phased according to the 
ephemeris of Spruit & Rutten (1998), whose binary phase 
zero is corrected to represent mid-eclipse of the white dwarf. 

The emission lines in VY Aqr vary considerably over the 
orbit. The H2O feature at 1.33/im also shows variation with 
phase. It does not appear to be present at phase <f> — 0.0. 
We have been unable to account for this variation; possible 
causes are discussed in section 5.1. There is no evidence for 



a change in the continuum slope in VY Aqr with orbital 
phase, with the slope always remaining within one standard 
error of the mean, given in section 4.2. 

The emission lines in WZ Sge show a strong variation 
throughout the orbit. There is a reduction in line strength 
at <f> « 0.9 and 4> « 0.0, corresponding to primary eclipse. 
The lines appear single peaked at <j) « 0.3 and <j> « 0.8. 

Time resolved spectra were not obtained for EF Peg 
and PX And. 

4.5 Skew mapping and Doppler tomography 

Skew maps were produced of both WZ Sge and VY Aqr. 
Skew mapping is a tomographic technique used to extract a 
value for the radial velocity of the secondary star, Kr. It is 
particularly useful when the spectral features are too weak 
for conventional cross-correlation techniques to work. It is 
described in detail in Smith, Cameron & Tucknott (1993) 
and has been successfully employed by Smith, Dhillon & 
Marsh (1998). We compared regions of the spectra unaf- 
fected by emission lines with the spectra of our M5 and M7 
dwarf stars. The first step was to remove the continuum from 
the spectrum of both the CV and the template stars. This 
is done by dividing by a first order polynomial fit, and then 
subtracting a higher order fit to the continuum. This ensures 
that line strength is preserved along the spectrum. Each 
of the time-resolved CV spectra were then cross-correlated 
with the template star spectra, giving a time-series of cross- 
correlation functions (CCFs) for each template. The skew 
maps were then generated by back-projecting the CCFs, 
in an identical fashion to time-resolved spectra in standard 
Doppler tomography (Marsh & Home 1988). If there is a de- 
tectable secondary star we would expect a peak at (0, Kr) 
in the skew map. 

Figure 4 shows the skew map of VY Aqr, made with 
G1866AB and a standard Doppler map of the Paschen-/? 
line in VY Aqr. When we first performed the skew mapping 
of VY Aqr, a strong peak appeared in the lower right-hand 
quadrant. Because we only cross-correlate with absorption 
features in the M-dwarf template, it is unlikely that the peak 
could be from anything but the secondary star, and so we 
believe that the unusual position of the peak is due to a dis- 
crepancy between the ephemeris of Thorstensen & Taylor 
(1997) and the actual ephemeris at the time of our obser- 
vations. Thorstensen & Taylor (1997) give their error in the 
period as ±0.00004 days, allowing for a discrepancy of up 
to seven orbits between the calculated and actual phases at 
the time of observation. Hence, both maps of VY Aqr have 
had a phase correction of -0.28 applied to the ephemeris of 
Thorstensen & Taylor (1997), in order to move the peak in 
the skew map to the expected position on the Kx=0 axis. 

The skew map of VY Aqr shows a significant peak, 
yielding a value of Kr = 320 ± 70 kms -1 . As a check, the 
time resolved spectra of VY Aqr were shifted to correct for 
the radial velocity of the secondary star and then averaged. 
The results are shown in figure 5. Secondary star absorption 
features are clearly visible in the corrected spectrum, which 
also shows a much stronger and sharper water band. We 
conclude that our value of Kr is appropriate, and not an 
artifact of the skew mapping process. Thorstensen & Taylor 
(1997) finds Kw — 49 ±4 kms -1 from Ha, giving a value for 
the mass ratio of the system of q = 0.15±0.04. This value of 
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Figure 3. Left-hand side: Time-resolved J-band spectra of the short period DN VY Aqr and spectra of the M-dwarf G1644C. Orbital 
phase is shown on the right-hand side of the plot, using the ephemeris of Thorstensen & Taylor (1997) — i.e without applying the phase 
offset described in section 4.5. Right-hand side: Time-resolved J-band spectra of the short period DN WZ Sge and spectra of the M-dwarf. 
Orbital phase is shown on the right-hand side of the plot, using the ephemeris of Spruit & Rutten (1998). In both plots the spectra have 
been normalized by dividing by the flux at 1.3 /Ltm. The spectra were then offset by a multiple of 1.0 for VY Aqr and 1.4 for WZ Sge. 
Also shown is the spectrum of an F8V star, normalized by dividing by the flux at 1.3 /im, which indicates the location of telluric features. 
The location of bad pixels is labelled by a bar showing the extent of the anomaly and the label 'bp'. 

q is consistent with the precessing-disc model for SU UMa should be noted that Kw, when derived from emission lines, 
stars (Whitehurst 1988), to which VY Aqr belongs. We are does not always reflect the true motion of the white dwarf, 
unable to derive component masses, as the inclination of 

the system is uncertain (Augusteijn 1994). Furthermore, it The Doppler map of VY Aqr shows a faint ring cor- 

responding to disc emission, and a peak we associate with 
the bright spot at (-500, -300) kms -1 . Assuming our phase 
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Figure 4. (Top) Left-hand panel: The Doppler map of the Paschen-/? line in VY Aqr. Right-hand panel: The skew map of VY Aqr, made 
with the M-dwarf G1866AB (M5V). The value of Kji is given by = + K y where Kx and Ky are found from the skew map by 
the position of the peak. Both maps have had a phase correction of -0.28 made to the ephemeris of Thorstensen & Taylor (1997). The 
predicted position of the secondary star, the path of the gas stream (lower curve) and the Keplerian velocity at the gas stream (upper 
curve) are marked on both maps, using the value of q derived in the text. The three crosses on the maps are, from top to bottom, the 
centre of mass of the secondary star, the system (at zero velocity) and the white dwarf. The circles plotted on the gas stream represent 
(from left to right) increasing distance from the white dwarf (in steps of one-tenth of the distance to the inner Lagrangian point, L\). 
(Bottom) Left-hand panel: The Doppler map of the Paschen-/? line in WZ Sge. Right-hand panel: The skew map of WZ Sge, made with 
the M-dwarf G1866AB (M5V). The predicted position of the secondary star, the path of the gas stream and the Keplerian velocity at 
the gas stream are marked on both maps, along with the centre of mass of the system components, in an identical fashion to that above, 
using the system parameters of Spruit & Rutten (1998). 



offset is correct, we find that the bright spot in VY Aqr 
is located downstream from the stream impact region. This 
is consistent with the optical and infra-red doppler maps 
of WZ Sge (see below) and with infra-red eclipse maps 
of IP Peg, which show that the infrared continuum bright 
spot can also appear downstream of it's optical counterpart 
(Froning et al. 1999). Work is in progress to determine if 
this is a common feature in infra-red accretion disc maps 
(Littlefair & Dhillon 1999). This consistency with other in- 



frared maps, as well as consistency with the radial velocity 
curves and light curves presented later in this paper gives 
us confidence that our phase correction is appropriate. 

Figure 4 shows a standard doppler map of the Paschen- 
P line in WZ Sge, along with a skew map of WZ Sge, made 
with G1866AB. The doppler map shows faint accretion disc 
structure and is dominated by emission from the bright spot 
at (-800, -200) kms -1 . Spruit & Rutten (1998) found that 
the optical bright spot was located downstream from the 
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+ CO - 



X} 
0) 



A. (/J.m) 

Figure 5. Average J-band spectra of the short period DNc 
VY Aqr and the M dwarfs G1644C and G1866AB. The upper 
spectrum of VY Aqr has been corrected for the radial velocity of 
the secondary star and represents an average spectrum in the rest 
frame of the secondary star. The spectra have been normalised 
by dividing by the flux at 1.3 fj,m and then offset by adding a 
multiple of 0.9 to each spectrum. Also shown is the spectrum of 
an F8V star, normalised by dividing by the flux at 1.3 /im, which 
indicates the location of telluric features. 



stream impact region. Using their ephemeris we find that 
the infra-red bright spot is not aligned with it's optical coun- 
terpart but is located further downstream from the stream 
impact region, consistent with our results for VY Aqr. The 
skew map of WZ Sge shows no peak from the secondary 
star but a bright streak corresponding to a large value of 
the cross-correlation function at (f> « 0.9. It is possible that 
this is due to correlation with secondary star features, es- 
pecially as it occurs near eclipse, although the absence of a 
strong correlation at other phases, and the fact that it is not 
exactly at <f> — makes this unlikely. 



4.6 Light Curves 

Figure 6 shows equivalent width (EW) light curves for 
WZ Sge and VY Aqr. The light curves have been folded over 
two complete orbits for clarity. They show a large amount 



of scatter, due to noise in the spectra, and hence should 
be interpreted with caution. Both the Hel (10830 A) and 
the Paschen-lines light curve in WZ Sge show an eclipse at 
0^0. 

The VY Aqr light curves show a systematic variation in 
brightness, with a minimum brightness occurring at <f> w 0.9 
in the Hel lines and at <f> £s 0.75 in the Paschen lines, where 
we have applied the correction of -0.28 to the ephemeris 
of Thorstensen & Taylor (1997). Such a variation could be 
produced either by obscuration of the bright spot by an op- 
tically thick disc, or by self-obscuration of an optically thick 
bright spot whose outer edge is brighter than the inner. 
Doppler maps of the Hel emission in VY Aqr show emis- 
sion from a bright spot which is not aligned with the Pa/3 
emission. This explains the lag between the Hel light curve 
and the Paschen-line light curve. The presence of any varia- 
tion at the orbital period in the light curves is surprising, as 
Patterson et al. (1993) found no variations on the orbital pe- 
riod in the optical continuum light curve of VY Aqr. Further 
infrared photometry seems desirable. 

4.7 Radial velocity curves 

The time resolved spectra were continuum subtracted and 
binned onto a constant velocity interval scale. In order 
to measure the radial velocities, we applied the double- 
Gaussian method of Schneider & Young (1980) to the 
Paschen-/3 line. This technique is sensitive mainly to the 
motion of the line wings and should therefore reflect the 
motion of the white dwarf with the highest reliability. The 
Gaussians were of width 500 kms" 1 (FWHM) and we varied 
their separation from 600 to 2400 kms -1 . We then fitted 



V — 7 — K sin( 



(1) 



to each set of measurements, omitting the points at primary 
eclipse. Examples of the radial velocity curves obtained for 
VY Aqr and WZ Sge are shown in fig. 7. 

The results of the radial velocity analysis are displayed 
in the form of a diagnostic diagram in figure 8. By plot- 
ting K, its fractional error ctk/K, 7 and <f>o as functions 
of the Gaussian separation it is possible to select the value 
of K which most closely matches Kw (Shatter, Szkody & 
Thorstensen 1986). If emission is disc dominated, one would 
expect K to asymptotically approach Kw when the Gaus- 
sian separation becomes sufficiently large. Furthermore, one 
would expect 4>o to approach 0. 

For WZ Sge it can be seen that the phase shift 0o does 
indeed fall towards with increasing Gaussian separation. 
The value of K steadily drops, reaching its lowest value at 
around K = 85 ± 20 kms" 1 . At this point the phase shift is 
0.17 ± 0.06 and not zero, as is typical for this sort of mea- 
surement in CVs (see e.g. section 2.7.6 in Warner 1995 for 
a discussion). This phase shift introduces doubt about the 
validity of adopting this figure for Kw- Several authors have 
attempted to measure Kw in WZ Sge with optical values 
ranging from 40-70 kms" 1 for Ha and 60-80 kms" 1 for H/3 - 
see Skidmore (1999). Spruit & Rutten (1998), using a light 
centres method which attempts to correct for contamina- 
tion by light from the hotspot, find a value for Ha of 40 ± 10 
kms" 1 . K-band infrared studies of WZ Sge yield a value for 
K, from the Br7 line, of 108 ± 15 kms" 1 (Skidmore 1999). 
Our value is compatible with the higher optical values. 
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Figure 6. Emission line equivalent width (EW) light curves of WZ Sge and VY Aqr. EWs were calculated from the Hel (10830 A) line 
and from the sum of the Paschen lines (P7 + P/3). The Hel line was separated from the Paschen-7 line by measuring the EW up to a 
point at the middle of the blend. Hence the Hel light curve is contaminated with a small amount of light from Paschen-7 an< l vice- versa. 
The light curves have been folded over two complete orbits for clarity. Phases for WZ Sge are from the ephemeris of Spruit & Rutten 
(1998), phases for VY Aqr are shown using the ephemeris of Thorstensen & Taylor (1997), with a correction of -0.28 applied. 




1 -1 
Orbital phase 

Figure 7. Radial velocity curves of Paschen-/3 for WZ Sge and VY Aqr, measured using a double-Gaussian fit with Gaussian separations 
of 1600 kms -1 (VY Aqr) and 2000 kms -1 (WZ Sge). Points at primary eclipse are omitted from the plot. The curves are folded over 
two complete orbits for clarity. Also shown are the fits to the curves (solid line) and the systemic velocities (horizontal dashed line). 
Phases for WZ Sge are from the ephemeris of Spruit & Rutten (1998), phases for VY Aqr are shown using the ephemeris of Thorstensen 
& Taylor (1997), with a correction of -0.28 applied. 



For VY Aqr the phase shift does not fall to zero, instead 
remaining roughly constant at 0o ~ 0.16, again typical for 
this type of measurement. The value of K remains fairly 
constant at K = 150 ± 50 kms -1 . This value is inconsistent 



with the value obtained by Thorstensen & Taylor (1997) of 
K w = 49 ± 4 kms" 1 for Ha. 

It is not easy to read reliable values of Kw from our 
diagnostic diagram. Thus, in an attempt to obtain reliable 
values for Kw for both VY Aqr and WZ Sge we employed 
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Figure 8. Diagnostic diagrams for WZ Sge and VY Aqr based on the double-Gaussian fits to Paschen-/3. Phases for WZ Sge are from 
the ephemeris of Spruit & Rutten (1998); phases for VY Aqr arc shown using the ephemeris of Thorstensen & Taylor (1997), with a 
correction of -0.28 applied. 



a modified version of the light centres method, described 
by Marsh (1988). In the co-rotating co-ordinate frame, the 
white dwarf has velocity (0,-Kw), and any symmetric emis- 
sion (e.g. from the disc) would be centred at that point. By 
plotting Kx = —K sin 4>q against Ky = —Kcoscjio for the 
different radial velocity fits, we find that the points move 
closer to the Ky axis with increasing Gaussian separation. 
This is because at large separations the fits are dominated 
by the line wings, which should be formed in the symmet- 
ric, inner disc. An extrapolation of the points to the Ky 



axis should give a measurement of Kw ■ Such a light centres 
diagram is presented in figure 9. 



The light centres in both plots seem to cluster at the 
edges of the bright spots as the gaussian seperation in- 
creases. Hence we are unable to extrapolate to Kx = and 
therefore unable to estimate Kw for WZ Sge and VY Aqr. 



© 0000 RAS, MNRAS 000, 000-000 



WZ Sge 



J -band spectra of CVs 11 
VY Aqr 
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Figure 9. Light centres of the Paschen-/3 emission in WZ Sge and VY Aqr, superimposed on the Paschen-/? Doppler maps. For clarity, 
strong emission is represented as dark areas for WZ Sge, and light areas for VY Aqr. Points are plotted for double Gaussian fits with 
separations ranging from 600 kms -1 (marked "low") to 2400 kms -1 (marked "high"). The dashed line represents Kx = 0, and is where 
the white dwarf is expected to lie. The centre of mass of the white dwarf (marked "o") is plotted using the values for Kw of Spruit & 
Rutten (1998) for WZ Sge and Thorstensen & Taylor (1997) for VY Aqr. Phases for WZ Sge are from the ephemeris of Spruit & Rutten 
(1998); phases for VY Aqr use the ephemeris of Thorstensen & Taylor (1997), with a correction of -0.28 applied. 



5 DISCUSSION 

5.1 Detection of secondary star features 

We consider here the reality of the secondary star detection 
in VY Aqr and EF Peg. The secondary star detection is 
based on three points of evidence; the continuum slope, the 
presence of the water band at 1.33 fim and the success of the 
skew mapping technique. We consider these points in turn. 

(i) VY Aqr has a red continuum. The continuum shows 
a change in slope around 1.1 /jm, a feature shared by the 
M-dwarf spectra. In fact, the continuum shape in VY Aqr 
strongly resembles that of the M-dwarfs. As the accretion 
disc and white dwarf are expected to have a blue continuum, 
it is reasonable to ascribe the red continuum seen in VY Aqr 
to light from the secondary star. This does not necessarily 
imply that the secondary dominates the IR light, as the 
amount of secondary star light needed to create an overall 
red continuum is strongly dependant upon the slope of the 
secondary star spectrum. The continuum slope in VY Aqr 
appears to show no variation with orbital phase, always be- 
ing within one standard error of the mean. 

(ii) The water band at 1.33 /jm lies in a region affected 
by telluric features (see fig. 1), so it is possible that incorrect 
telluric correction could affect this feature. In order to try 
and estimate the extent of this effect, spectra of VY Aqr 
were calibrated using each of the F-stars observed. It was 
seen that there was excess absorption at 1.33 /jm, even 
when the stronger telluric features around 1.14 /jm were 
over corrected. Therefore, incorrect telluric correction can- 



not account for the observed feature. The presence of the 
unexplained variations in strength of the water feature with 
orbital phase (as shown in figure 3) is disturbing, however. 
Following the method outlined above, we find that incorrect 
telluric correction cannot account for this variation. Also, by 
repeating the reduction process with a range of parameters 
and methods, we determine that systematic errors in the 
reduction process do not account for the variation either. 
The feature is only present in the time resolved spectra at a 
significance of one standard error, so we might attribute the 
variation in the band to random noise. However, the detec- 
tion of this feature cannot be regarded as certain. It is noted 
here that the percentage contributions of the secondary to 
the J-band light in VY Aqr and EF Peg are dependant on 
this feature, and are hence also uncertain. 

(iii) The skew mapping technique correlates with a sec- 
ondary star template. Strong cross-correlations are obtained 
with this method. When the time-resolved spectra are aver- 
aged in the rest frame of the secondary star, absorption fea- 
tures are stronger and sharper. Further, the skew mapping 
suggests a correction to the ephemeris for VY Aqr, which is 
consistent with the radial velocity curves and light curves. 
Thus we are confident that the skew mapping is correlating 
with real secondary star features. 

None of these items, considered alone, provide concrete 
signs of the secondary star. However, in combination they 
provide a strong body of evidence. Hence we are confident 
that the secondary star has been detected in VY Aqr. The 
detection in EF Peg is based only on the continuum shape 
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and the water feature at 1.33 /im. Hence the detection of the 
secondary in this case is less reliable than that in VY Aqr. 



6 CONCLUSIONS 

(i) The relative ease of observation in the J-band, and 
strong absorption features of late-type dwarfs in this spectral 
region make J-band studies highly desirable for CVs with 
undetected secondary stars. 

(ii) The secondary star has been detected in VY Aqr and 
EF Peg, although the detection was too weak, and insuffi- 
cient field-dwarf templates were observed, to make an esti- 
mate of the spectral type. 

(iii) The contribution of the secondary star to the total 
J-band flux was found to be 65% for EF Peg and 70% for 
VY Aqr. Upper limits of 10% and 20% were found for the 
contribution of the secondary to the total J-band flux for 
WZ Sge and PX And respectively. 

(iv) A value of Kb, = 320 ± 70 kms -1 was found for 
VY Aqr and the mass ratio was calculated as q = 0.15±0.04, 
assuming Kw, as derived by Thorstensen & Taylor (1997) 
is correct. 

(v) The spectral type of WZ Sge was constrained to be 
later than M7.5V. 
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Table 1. Journal of observations. Each spectrum consists of 240 s total exposure time. Gaps between objects were used for arcs, fiats 
and F-star spectra. An object is listed twice where such a gap broke continuous observation. The long gap between the observation 
of G1866AB and WZ Sge is due to cloud. Observations of VY Aqr between 07:20 and 08:54 contained no useful information, due to 
very poor seeing. Spectral types are those given by Henry, Kirkpatrick & Simons (1994). Orbital phase is shown for WZ Sge using the 
cphcmcris of Spruit & Ruttcn (1998), PX And using the ephemeris of Still, Dhillon & Jones (1995) and VY Aqr using the ephcmeris of 
Thorstensen & Taylor (1997). The orbital period for EF Peg is estimated from the supcrhump period (Howell et al. 1993). 



Object 


•^central 


Date 


UT 


UT 


No. of 


Period 


Phase 


Phase 


Spectral 


J-band 




(H 




start 


end 


Spectra 


(days) 


start 


end 


Type 


mag. 


WZ Sge 


1.175 


08/08/98 


06:40 


08:33 


24 


0.0567 


0.24 


1.43 




15.3 


WZ Sge 


1.175 


08/08/98 


08:47 


10:36 


24 


0.0567 


1.80 


2.97 




15.3 


G1866AB 


1.175 


08/08/98 


10:56 
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4 








M5V 
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14 


0.0567 
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14.2 


G1644C 
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06:20 
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08:54 


20 


0.0635 


0.81 


1.19 




14.2 


VY Aqr 


1.175 


09/08/98 


09:22 


10:53 


20 


0.0635 


1.50 


2.36 




14.2 


EF Peg 


1.175 


09/08/98 


11:36 


13:09 


20 


0.0837 










EF Peg 


1.175 


09/08/98 


13:33 


14:46 


16 


0.0837 
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Table 2. Wavelengths, equivalent widths and velocity widths of the most prominent lines in the J-band spectra of the novalike variable 
PX And, the DNe WZ Sge, EF Peg, VY Aqr and the M-dwarfs G1644C and G1866AB. The line identifications were made using the lists 
and spectra of Kirkpatrick et al. 1993 (marked a ), Joyce et al. 1998 (marked b ), Jones et al. 1994 (marked c ), Jones ct al. 1996 (marked 
d ), Tinney, Mould & Rcid 1993 (marked e ), Lang 1986 (marked i) and the UKIRT online documentation (marked 9 ). Equivalent widths 
for the water band around 1.34 fim were measured in the range 1.32—1.34 (im. The vertical bars following lines of similar wavelength 
indicate that the measurements apply to the entire blend. Wavelengths for VO, TiO, H2O and FeH refer to the band heads, or the 
edge of the band, if headless. The two-letter codes indicate that the line is either not present (np) or that the line is present but is not 
measurable (nm). 
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Table 3. The Contributions of the secondary star to the total J-band flux. 



Object Secondary Error M-dwarf Mag. of 

Contribution used Secondary 



WZ Sge 


<10% 


5% 


M7 


>16.2 


PX And 


<20% 


5% 


M5 


>15.4 


EF Peg 


65% 


15% 


M7 


16.7±0.5 


VY Aqr 


70% 


15% 


M5 


15.7±0.5 
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